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Conformations of isolated homo- dendrimers of G — 1 — 7 generations with D — 1 — 6 spacers 
have been studied in the good and poor solvents, as well as across the coil-to-globule transition, by 
means of a version of the Gaussian self-consistent (CSC) method and Monte Carlo (MC) simulation 
in continuous space based on the same coarse-grained model. The latter includes harmonic springs 
between connected monomers and the pair-wise Lennard-Jones potential with a hard core repulsion. 
The scaling law for the dendrimer size, the degrees of bond stretching and steric congestion, as well as 
the radial density, static structure factor, and asphericity have been analysed. It is also confirmed 
that while smaller dendrimers have a dense core, larger ones develop a hollow domain at some 
separation from the centre. 



PACS numbers: 36.20.-r, 36.20.Ey, 61.25.Hq 



I. INTRODUCTION 

Dendrimers represent a class of well defined hyper-branched macromoleculestLwhich can be synthesised vi, 
sequence of carefully controlled repetitive reactions producing regular structures! There are high aspirations! 
of chemists and piwsicists regarding applications of such novel polymers as advanced materials for medicinalB uses, 
superior catalystsu'l-1, and drug delivery vchiclesa, to mention just a few. Importantly for such applications, dendrimers 
are quite flexible molecules, which have a well accessible interior part as well as a large exterior surface. Furthermore, 
dendrimers can controllably change their size and density distribution under the influence of external conditions. Some 
of the potential .applications of dendrimers moreover require a hollow interior, which could be used for accommodating 
guest moleculesu. In addition, spacers, charges, and distinct types of monomer units should allow one to further fine 
tune particular behaviour of dendrimers. 

Although the ideal model of dendrimers with the harmonic springs or random walks can be solved analytically 
such results are of a limited value since the excluded volume interactions are of paramount importance for dendrimer 
However, including these in order to describe the range of the coil-to-globule transition in sufficient detail is a 
considerable challenge even for the homo- dendrimer. This problem has attracted a number of researchers yielding a 
significant body of theoretical and computational work in recent years, with some experimental resultsc£l available to 
assist them. _ 

De Gennes and HervetEj have performed a seminal calculation using a version of the Edwards' self-consistent field 
method upon an assumption that an onion-like layered distribution of the dendrimer generations applies. They 
have predicted a density minimum at the core of the dendrimer and the scaling law for the radius of gyration, 
TZ g ~ TV 1 / 5 , whose N is the total number of monomers, both of which henceforth have been debated. Later on Boris 
and Rubinsteinli-3 have proposed a more complicated self-consistent mean-field and related Flory-type theories. They 
have found that the density is monotonically decreasing with a maximum at the core. However, such a theory does 
not contain spatial correlations and thus may have considerable limitations in predicting the internal structure of 
dendrimers. 



* Corresponding author. Web: http://darkstar.ucd.ie; E-mail: Edward. Timoshenko@ucd.ie 



1 



A version of the Gaussian variational theoryQ, which relies on a somewhat problematical virial representation of 



the excluded volume interactionsEOllj, 
viscosity of dendrimers in a good solvent 



s been used by Ganazzoli et al for describing the size, shape and intrinsic 
Based on this approach-the dynamic properties such as the relaxation times. 



spectrum, dynamic structure factor and the viscoelastic modulij£3, as well as the inter-chain theta-temperaturecll 
have been also investigated. Unfortunately, constraints of the numerical procedure employed limited such studies to 
dendrimers with D = 1, 2 spacers only. 

A number of simulation studies using either Molecular Dynamics (MD) or Monte Carlo (Mil), both on and off- 
lattice, whether atomistically detailed or coarse-grained, have been also performed. Naylor et a£3 have discussed the 
MD simulation of PAMAM (poly(amido-amine)) dendrimers. This study, albeit generally interesting, could cover a 
timescale of hundreds of picoseconds only. Later Lescanec and Muthukumaio have proposed a kinetic growth model 
in a 3-d off-lattice MC. They have reported the highest density at the core and have also found that the terminal 
monomers can traverse the dendrimer. Their approach, however, was in essence a non-equilibrium one as sufficient 
equilibration was not been allowed. This may explain a rather small reported value of the swelling exponent v ~ 0.22 
in the proposed scaling law lZ g ~ N v D 1 / 2 for the good solvent. Next, Mansfield and Klushinc-3 have performed a MC 
simulation on a diamond lattice with D = 7 spacers modelled as self-avoiding walks. Their paper has managed to 
overcome some of the previous limitations and it contained some important insights. Nevertheless, the D dependence 
has not been investigated there. Furthermore, a possibility of lattice artefacts for.-the restricted branch points may 
raise some doubt in the more refined features of their predictions. Chen and CuiEa have performed off-lattice MC 
simulations for a good solvent in a model with freely rotating bonds. They have also attempted to fit the radius of 
gyration as R g ~ D" 2^ 2l/ -^ G G 1 "" , i.e. in terms of both D and the number of generations G. While the value of v has 
been related to the inverse fractal dimension of the open/ring chain, just as in Ref. [2^ containing the renormalisation 
group calculation by Biswas and Cherayil, the meaning of the G-dependent part of the scaling law lacked such clarity. 

Further MD simulation by Murat and GrestEO has relied on a coarse-grained model with the Lennard-Jones non- 
bonded and the FENE bonded interactions employing a white noise for the solvent and temperature effects. This 
simulation involved sufficient equilibration times and has yielded a number of interesting results. For instance, the 
dendrimer size at a fixed number of spacers D in a good solvent was found to scale as R g ~ N v with the value of the 
exponent v surprisingly close to 1/3. Interestingly, the indication about a hollow domain inside the dendrimer, which 
has been first noted in Ref. |4| has been reiterated and further elaborated by performing the partial densities analysis 
over generations. These authors have also looked at the dendrimer contraction on changing the solvent quality from 
good to poor. They have concluded that the law lZ g ~ N 1 ^ 3 (at fixed D) applies for the more compact globular 
dendrimer, which seemed reasonable from the point of view of the space filling argument. 

Despite this considerable progress, understanding of the homo- dendrimers is not as yet fully satisfactory. First, a 
reasonably large span of spacer and generation numbers under variable solvent conditions has until now been hard to 
study. Second, it was difficult to achieve both a reliable equilibration and a good averaging statistics afterwards at the 
level of affordable computational times. Clearly though, these requirements must be met at the same time in order to 
be able to extract reliable scaling predictions and to make firm conclusions on the intricacies of the internal structure 
of dendrimers. Moreover, not only the good, but the theta- and poor solvents, as well as the thermodynamics of the 
coil-to-globule, are of interest. Perhaps, particular deficiencies of the techniques employed, an insufficient accuracy, 
and a limited range of systems may explain some of the discrepancies which remained between the results of different 
approaches. Therefore, we believe that it would be important to revisit the problem. 

Here we shall try to re-examine the fundamentals of dendrimers in a broader context of their overall conformational 
behaviour as a function of the solvent quality. On the one hand, this may help to resolve some of the remaining 
controversy in the concrete predictions. On the other hand, it may shed some new light on the problem thanks to 
studying the observables which hitherto have not been investigated and to performing a type of analysis which is 
more novel. Last but not least, one-piay hope to improve the statistics and to study a broader range of dendrimers 
by using a well optimised MC codeEao on faster modern computers^. 

Our recent progressed on extending the Gaussian self-consistent (GSC) method may offer some advantages and 
permit comparisons to be made with the simulation data. The improved technique avoids the use of a virial expansion, 
but instead operates directly with a given Hamiltonian involving any pair-wise bonded and non-bonded interactions, 
just as in direct simulations. Ref. [32| deals with all the technical details and it contains a comprehensive comparison 
with the MC data performed for open, ring and star homopolymers with variable flexibility in a good and poor 
solvents. The method appears to be relatively accurate for a number of observables and its limitations are well 
understood. Because the new GSC technique is computationally very fast and since it directly yields the equilibrium 
averages for the observables of interest, it will allow us to consider a larger range of system sizes and a broader region 
of solvent conditions than with the equivalent MC simulation. 
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II. MODEL AND NOTATIONS 



The current coarse-grained homo- dendrimcr model is based on the following Hamiltonian (energy functional )0@ E1 
in terms of the monomer coordinates, X^: 

^^EM^-X^ + i £ [/(WflX.-X.I). (1) 

Here the first term represents the connectivity of the dendrimer with harmonic springs of strength Kij introduced 
between any pair of connected monomers (which is denoted by i ~ j). The second term represents pair-wise non- 
bonded, specifically van der Waals interactions between monomers. We shall adopt the Lennard- Jones form of the 
potential, 

!+oo, r < 

with hard core part and the monomer diameter S°\ where is the dimensionless strength of the interaction. 

It is important to introduce convenient notations. Let us consider a dendrimer of G generations, each of which 
consists of D spacers, with F functional core (corresponding to the generation g = — 1) and F functional other branch 
points (see Fig. [I]). Thus, if D = a and G — b we shall briefly denote such a dendrimer as DaGb as in Refs. [IJ]|20[ 
Then the total number of monomers will be, 

(F-l) G+1 - 1 

N=l + F D [ - . (3 ) 

Note that below we shall consider results for a particular example of tri-functional branching F = Fq = 3. It 
is convenient to represent monomer indices via triads i — (g,d,tp g ), where g = 0,1,... G is the generation index, 
d = 0, 1, . . .D — 1 is the spacer index, and ip g = 0, 1, . . . Fq(F — I) 9 — 1 is the angular in-generation index. It is also 
useful to define the shell radial index via p = Dg + d + 1, and hence instead of the triad, one can use the 'polar 
coordinates' index notation i = (p,ip g ). Both of these will be used interchangeably henceforth. 



III. TECHNIQUES 
A. Monte Carlo technique 

We use the Monte Carlo (MC) technique with the standard Metropolis algorithmEl and local monomer moves, based 
upon the implementation described by us in Refs. p8|]3C| . The position of a randomly picked monomer in the spherical 
coordinates is displaced by Ar = 5r\, A8 = n0i, A<p = 2tt4>i where (r*i,0i,<^i) is a triple of independent standard 
uniform deviates and 5 is an additional parameter of the MC scheme characterising the timescale involved in the Monte 
Carlo sweep (MCS), the latter being defined as TV attempted MC steps. The Metropolis check, AE < — fcsThiri, 
where r± is yet another standard uniform deviate, is then used as a criterion for accepting trial new conformations. 
As a result, the ensemble of independent initial conformations would converge to the Gibbs equilibrium distribution 
with increasing number of sweeps. Initial conformations of dendrimers were taken as planar and rather overstretched 
ones (akin to the schematic representation in Fig. Q), but then they were subjected to extensive equilibration for a 
required time before any simulation was commenced. 

To ensure good equilibration, the behaviour of global observables such as the energy and radius of gyration was 
monitored. Upon reaching equilibrium these cease to have a global drift and start exhibiting characteristic fluctuating 
behaviour around well defined mean values. Then, depending on the dendrimer size, about Q = 2 • 10 4 of statistical 
measurements have been performed, typically separated by about 40 TV 2 of attempted Monte Carlo steps to ensure 
statistical independence of sampling. The mean value and error of sampling of an observable A are then given by the 
mean (A) = (1/Q) A and by ±^J{AA) 2 /Q respectively. 
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B. The Gaussian Self— Consistent method 



The main objects in the GSC method are the mean-squared distances between monomers, 

©«(<) = ^(XiW-X^t)) 3 ^. (4) 

The GSC method is based on replacing the stochastic ensemble for Xi with the exact Hamiltonian in the Langevin 
equation of motion onto the trial ensemble xf^(t) with a trial Hamiltonian H^(t), which is a generic quadratic form 
with matrix coefficients called the time-dependent effective potentials, 

fr(°)[x(t)] = i^;^(t)x < (t)x i (i). (5) 

y 

Then one requires that the inter-monomer correlations satisfy the condition, 

(X i (t)X j (t)) = (x[°\t)xf\t)) 0) (6) 

which means that the trial ensemble well approximates the exact one. This yields the self-consistent equationscl, 
which in the absence of the hydrodynamic interactions are, 

^V ij{ t) = - 2 - VWt) - V jk (t)) ( ^« - MWM) . (7) 

Here is the friction coefficient of a monomer and the instantaneous free energy has the same functional expression 
via the instantaneous T>ij (t) as it has at equilibrium. 

The expression for the free energy, A = £( tot ) — TS has been given in Ref. |3^. We would not like to reproduce 
that expression here as it is somewhat cumbersome. However, it should be mentioned that the total mean energy 
includes the energy of the bonded and non-bonded interactions £ as well as the hard-core, contribution as follows, 
£(tot) = £ + j\i hs \ The second term is given by the generalised Carnahan-Starling terrr£3i£3, 

k B T -2^ {\~ ni y ' 7,1 -L* d (o) - W 




F W[y} = V^JIZ (9 ) 



where r\i is known as the packing coefficient of the monomer i, which would have an important role below. 

The stationary limit of the GSC Eqs. (Q) produces the equations for the minimum of the free energy, which are the 
same as those derived from the Gibbs-Bogoliubov variational principle. Although here we shall only be concerned 
with the equilibrium, the numerical solution of Eq. (Q) , applied until the stationary limit is reached, presents one of 
the most efficient technioues for finding the global free energy minimum. This, based on the fifth order adaptive step 
Runge-Kutta integratoilij, is used for obtaining the results from the GSC technique. 

It is also important to map out all equivalent pairs of the matrix 2?y , which is discussed in Appendix ^ based on 

the concept of the topological tree, so that only C^ d = (G + I )(D /2){D(G 2 /3 + (13/6)G + 2) + 2 + G/2) independent 
elements out of total Ctot — N(N — l)/2 remain in it due to the kinematic symmetries of the dendrimer. It is 
interesting to remark, more generally, that exactly the same symmetries will exist for a hetero- dendrimer as long as 
the monomers at each shell index p are identical to each other. 

The kinematic symmetries for the matrix T>ij also yield analogous symmetries for the matrix Vij. Thus, the 

computational expenses per step in our calculations are of order t c ~ N G^, where cf^ d ~ D 2 G 3 /6 is the total number 
of independent elements in the T>ij matrix (see Eq. (|Al[) ). These symmetries significantly reduce the computational 
times compared to the equivalent system in the MC method, where such symmetries only appear in the observables 
after averaging over the statistical ensemble. For comparison, the computational expenses per step in MC are of order 
t c ~ N At S, where At ~ N 2 is the number of MC steps needed to ensure good statistical independence between 
measurements, and S is the number of measurements needed for sampling of observables. Typical values of S should 
be of order of 10 4 — 10 5 for good accuracyE2l. For example, for D3G5_jdendrimer (consisting of N — 568 units) in 
the good solvent, — 1, the GSC method takes about 22 minutesU to reach equilibrium and produce its data, 
whereas the MC manages to make only 100 measurements of observables per hour after required equilibration. Thus, 
the GSC is about 10 3 times faster than the MC for producing the same results here. 
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IV. RESULTS 



We shall restrict ourselves to springs of equal strength, «y = 1, and choose the hard sphere diameter equal to the 
statistical length, — £, in Eq. (Q) as in Ref. ||(| Moreover, henceforth we shall use the mean energy £ expressed 
in units of ksT and the mean-squared distances T>ij and mean-squared radius of gyration, 31Zg, expressed in units 
of^ 2 . 

In our figures below we shall present data from MC simulation via thick lines and empty circles, whereas data 
from GSC theory via thin lines and small filled circles. In this work we have analysed the following dendrimer sizes: 
G = 1 - 7 for D = 1 - 3 and G = 1 - 6 for D = 4 - 6. 



A. Good solvent 



We shall start by examining the dendrimer overall size and structure in the good athermal solvent, C/^ ^ = 0. 
The mean-squared radii of gyration 37?.g of dendrimers with varying number of spacers D vs the number of 
generations G are depicted in Fig. ||. Clearly, 37?. 2 increases with both D and G, but in a rather different manner, 



which will be discussed in full detail below in Subsec. [VC. Corresponding MC and GSC curves follow each other 



quite closely, but the GSC theory increasingly overestimates the dendrimer size with increasing D and G, in a manner 
similar to our discussion in Ref. [32|. 

As G increases, dendrimers become more spherical, which is illustrated via the average normalised 'axes of inertia' 
\( a ) defined by Eqs. (9,10) in Ref. These converge towards a perfect sphere limit, A( a ) = 1/3, as G grows in Fig. ||. 
Interestingly, these quantities seem to be practically independent of the number of spacers D, which is illustrated by 
using circles (D = 4) and lines (D — 2) in the figure. 

Typical snapshots of dendrimers D3G3 and D3G7 in Figs. [| from MC simulation do confirm the observation that 
dendrimers become more spherical as G increases. Generally, a fairly large dendrimer tends to fill in the available 
space quite effectively and relatively uniformly, although some local dilutions and accumulations remain. The size of 
the dendrimer is well defined with very few branches outstretching beyond it. 

The shape of the static structure factor (SSF) in the Kratky representation in Fig. || for small G is reminiscent 
of that of a ring, or a star with few arms (see Figs. 9,10 in Ref. [30]), but tends more to that of a globule with 
increasing number of generations. Even though the tail of the plot still increases with p, it develops a characteristic 
oscillating behaviour. This is consistent with the trends in \( a > above, particularly so since the Independence is 
almost unnoticeable in S as well. 

The total and partial densities, 

/>oo 

gf)(r)=y^{5{vi-T cm -v)), 4tt / r 2 drg^(r)=N g = F D(F-iy. (10) 

of the dendrimers in Figs. |^ allow one to understand the monomer distributions in more detail. In a relatively small-G 
dendrimer (Fig. a) the total density monotonically decreases with the separation from the centre-of-mass r cm . The 
highest density at the origin also occurs for the partial density of inner generations, and it follows a similar pattern. 
Remarkably, the terminal generation density is rather delocalised spreading evenly from the outer boundary towards 
the very centre of the dendrimer. This can only be rationalised by allowing a possibility of the ends re-entering 
the dendrimer and reaching far inside, although the density of the terminal generation exceeds that of the others 
somewhat near the very edge of the dendrimer. 

These features are further complicated for dendrimers of a greater number of generations, typically for G > 5, 
at least in the considered cases of D = 2,4. While the behaviour of the terminal generation density remains very 
similar in Fig ^b, the density of non-terminal generations, and hence the total density, develop a well noticeable dip 
at certain separation. This indicates that larger-G dendrimers have the lowest density domain fairly close to the 
centre, but not at the centre itself, which always has the highest density of all. We may note here also that while the 
core monomer coordinate does differ from that of the centre-of-mass somewhat, they are quite close to each other 
in a large-G dendrimer. Thus, this most hollow domain is not located where the core monomer is and its separation 
from the core steadily increases with G. Clearly, from Fig |(jb one can see that the increase of the total density at the 
origin comes entirely from the density of the innermost generations. However, the core domain is at its most dense 
in a small-G dendrimer and it becomes less dense for larger dendrimers (in terms of both G and D actually) as can 
be seen by comparing the values of g^(r = 0) in Figs, [ga and b. This 'entropic pull' effect is due to an increased 
entropy thanks to extra branches attached to a larger dendrimer, leading to a higher stretching of the springs near 
the core with larger G. 
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Indeed, in Figs. |7]||[)]we depict the mean-squared distances between nearest neighbours along a branch vs the shell 
index p. This is perhaps the best way to characterise such stretching. While the main bodies of Figs. present the 
data from MC simulation, the insets exhibit the same results obtained from the GSC theory. Cleanly, the concrete 
numbers for P( p _ lj o) ( Pj o) are somewhat overestimated by the GSC theory compared to the MC dataE3, especially for 
small p in relatively large-G dendrimers. Nevertheless, the particular patterns and overall shapes of the MC and GSC 
plots are rather similar. 

As can be seen from Fig. [j], for small dendrimers D2G1-3, the spring from the core to the first monomer is stretched 
most, from the first to the second one is less, and so on, with the function 2?( p _i ! o) (p.o) decreasing monotonically. 
However, for dendrimers with G > 3 the function develops a characteristic descending step-like behaviour: the spring 
stretching within the same generation is about the same, but it has a dramatic drop when passing the branch point 
to the next generation. In all cases though, the bonds coming from the core are stretched most, with this feature 
becoming stronger for larger G, whereas the terminal bonds are the least stretched of all, with their length being nearly 
G- independent. The step-like behaviour can be easily interpreted by remarking that the number of springs doubles 
at any branch point. Hence the total unfavourable energy change due to a given magnitude of bonds stretching would 
be more significant near the edges than closer to the core. Note also that the step like dependence vanishes for the 
terminal generation. 

Plots in Fig. || for D = 4 spacers show that the stretching near the core increases significantly with G here as well. 
One can also see that the intra-generational steps develop a non-monotonic behaviour. Bonds within a non-terminal 
generation first become shorter with increasing index d until about D/2, then grow longer again when approaching 
the next branch point. This U-pattern occurs because there is more steric congestion near the branch points than 
in the middle of generations. The terminal generation is the least stretched here as for D = 2 and it has a similar 
monotonically diving behaviour ending up at a G independent stretching value. Similar plots for dendrimers with 
D = 6 spacers in Fig. |^ obtained from the GSC theory show the U-pattern becomes more pronounced as D increases. 

Finally, to address the issue of steric congestion we shall analyse the behaviour of the monomer packing coefficient 
i]i defined by Eqs. (^,^|). This quantity shows the volume fraction occupied by all other monomers around a given 
monomer number i, with the limit r\i — > 1 corresponding to the fully packed situation according to Eq. (|J). 

In Figs. [^,[n] the packing coefficients vs the shell index p along a branch are shown for different dendrimers with 
G = 6 and D = 3 respectively. For the former case it is convenient to normalise the horizontal axis as p/D so that 
branch points occur at integer values 0, . . . , G with G + 1 corresponding to the terminus. In Fig. [l^ different curves 
correspond to different numbers of spacers D. The curve for D = 1 gives a general outline of r\ p behaviour. The 
packing coefficient, and hence congestion, first increases, reaches a maximum and then rapidly decreases towards 
the termini, which are less crowded than the core domain. A larger number of spacers D leads to an U-pattern 
behaviour for g < G and a rapid dive for the terminal generation, similar to our discussion of Fig. ^. Again, the 
middle of generations, where rj p has smooth minima, has a lower steric congestion since the monomers have only two 
nearest neighbours along the chain there compared to three for branch points, where r\ p has sharp maxima. As for 
the terminus, it has a single nearest chain neighbour only, resulting in a pronounced minimum of rj p . 

In Fig. ^ we likewise present the packing coefficients for dendrimers with increasing number of generations G and 
the fixed number of spacers D = 3. The general behaviour of r\ p is similar here. For small G < 5 the function r\ p 
oscillates around a constant, dropping down dramatically on the terminal generation. However, in larger dendrimers, 
G > 5, the function r/ p starts from much lower values for small p and then increases overall up to penultimate 
generations, dropping down for the terminal g. Such behaviour is consistent with a lower density domain near the 
core seen via the spatial density distribution in Fig. ^b. For the D3G7 dendrimer, the lowest packing coefficient value 
occurs at around the middle of the g — generation (monomers 1,2,3 in diagram Fig. |l], which are buried deep inside 
the dendrimer as in Fig. |^b). Note also that here we have ?7_d(g+i) > Vo, indicating the onset of overcrowding in the 
terminal generation, which would become very severe for larger G. Overall, the packing coefficient typically changes 
nearly twice in its value for relatively large D and G = 6, whereas the range of change in r\ p increases further with 
increasing G. 

Limited simulations with a smaller monomer diameter d^> = 0.3 1 have been also performed and similar results 
have been obtained. Dendrimers with such reduced repulsions are more compact and, clearly, as is decreased 
further they change in a cross-over manner towards the behaviour of the ideal dendrimers. 

B. Coil— to— globule transition and poor solvent 

Now we shall bring our attention to the coil-to-globule transition. Upon increasing Lennard- Jones attraction C/' ' 
the dendrimer contracts and the mean-squared radius of gyration decreases as shown in Fig. |l2|a. The magnitude of 
this change is quite significant, being over 20 times for D6G5 dendrimer. As our plot is drawn in a single logarithmic 



G 



scale one can see that the region of maximal relative change in value of 31Z 2 , occurs at around ~ 2.5, with the 
change being more pronounced in larger dendrimers. Clearly, from this plot and from Fig. |l2|b, which depicts the 
derivative of the specific mean energy, one can conclude that the transition is continuous (second-order like). 

The quantity d(£/N)/dU^ ' changes from one asymptotic behaviour at small to another at large arguments, 
reaching eventually a constant value corresponding to the maximal packing. The maximal change in this quantity 

occurs at UgP ~ 2.5, which may be viewed as a reasonable practical definition of the transition point!!. As the 
transition is continuous it has a finite width and hence any such definition is somewhat arbitrary. However, as N 
increases the transition narrows and becomes sharper. It would be easy to consider the point of vanishing of the non- 
bonded part of the second virial coefficient (binary integral), u^j = (1/2) J dr(l — exp(— which occurs at 
about Uj^j ~ 1.386. Unfortunately, such a definition of the theta-point is unacceptable as it produces a result which 
does not depend on the polymer connectivity or size. Instead, one has to consider the total iv- 2 ' with the full account 
of the bonded and non-bonded interactions. This, however, is technically quite complicated as the Mayer-Ursell 
expansion is not easy to apply to dendrimers with both types of interactions involved. Thus, we shall use U £ , as a 
convenient numerical definition of a point which becomes close to the exact theta-point for sufficiently large N . 

In Tab. | we report the values of Ujh for dendrimers with varied D and G. Clearly, either increased number of 
generations or spacers, and hence increased N, would lead to a lower value of U £ , (larger temperature T £ >) of the 
transition. Indeed, the point U^) corresponds to the compensation of the positive conformational entropy S^ gau ^ 

with the negative two-body £^ energy growing as iV 2 and much faster than S^ gau ^ in N, so that smaller U £ °) would 
suffice for the compensation. More non-trivially though, if we compare, for example, D3G6 and D6G5 dendrimers, 
both having about the same N, the larger D system would have a lower U £ , (larger temperature Tg')- 

For linear or ring polymers one can alternatively define the theta-point of the current model by finding a value 
of at which the radius of gyration of the system matches that of the corresponding ideal system (g^ ' =0). In 
Tab. such values are presented in parentheses. Clearly though, this definition is inappropriate for dendrimers 
as u!£ shifts deep into the globular area for dendrimers with sufficiently large G. 

The region beyond Ug/ corresponds to a globular conformation of the dendrimers, to discussion of which we shall 
now turn. In Fig. [l^ we display the mean-squared radii of gyration for dendrimers with different number of spacers 
D vs the number of generations G. The solid lines correspond to the globule with [A ) = 7, that is, well beyond the 
transition region, whereas the dashed lines correspond to similar ideal dendrimers (d^ = and U^ = 0). These 
grow with D at a fixed G, and vice versa in both cases. Interestingly, while the globule is more compact than the 
corresponding ideal dendrimer at relatively small G, for D1G>6 this trend reverses. This is due to a significant role 
of the repulsive part of the non-bonded potential in dendrimers with excluded volume interactions and it shows a 
limited relevance of the model of ideal dendrimers when applied to the coil-to-globule transition. In fact, the scaling 
behaviour at the transition point differs from that of the ideal dendrimer unlike in the simpler cases of open or ring 



polymers. Such issues will be examined by us further in Subsec. [V C 



Next, we shall compare the mean-squared distances between nearest neighbours along a branch vs the shell index 
p for the globule in Fig. |l4| to that in the good solvent in Fig. ^. The function 2?( p -i,o) (p,o) here has a similar 
U-pattern inside each generation, a dive at the terminal generation, and it increases as G is increased, just as before. 
However, the descending stairway behaviour seen for the good solvent is no longer present now even for large G, and 
the function oscillates around a constant value. Thus, the bond stretching near the core is reduced, and apart from 
the intra-generational variation and the special role of the termini, the globular dendrimers are fairly evenly stretched 
as well as have uniform density (see inset in Fig. ||a). Note also that the maxima occur at pairs of subsequent chain 
indices p. This means that the mean-squared distances from the pre-branch monomer to the branch point are almost 
the same as from the branch point to the post-branch one, which indicates similarity in overall bond stretching with 
increasing p and hence an approximate equivalence of all branch points in the globule. Therefore, the connectivity 
structure of the dendrimer is less significant and more locally manifested in the globular state. 

Similar conclusions can be made also from the behaviour of the monomer packing coefficient for globular dendrimers 
in Fig. [15]. Again, rj p only depends periodically on the spacer index d, being nearly identical for all generations but 
the terminal one. The peaks occur at the branch points and the termini have the lowest packing coefficient of all just 
as for the good solvent in Fig. |ll|. Note, however, that the actual magnitudes of i] p are much larger for the globule 
than for the good solvent and approach 1/2. Moreover, the mean value of r] p in the globule steadily increases with 
N and there is no hollow domain in the globular conformation even for G > 5. Indeed, the total density g^'(r) of a 
globule dendrimer shown in the insets of Figs. ^a,b is practically uniform and rather localised within the size of the 
globule, thus having higher values than for the good solvent. Remarkably, the terminal generations still can traverse 
the dendrimer even in such a dense state. 
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C. Scaling relations 



First, one can fit the mean-squared radii of gyration at a fixed number of generations G in terms of N expressed 
via D by Eq. (0) as, 



Good solvent: A /3ft2 = a G N VG , (11) 



Theta-solvent: J3K 2 = a G N VG , (12) 



Poor solvent: ^3K 2 g = A G N VG . (13) 

The values of the exponent vg obtained from fitting of the data are presented in Tab. H, 

MC: v G ~ 0.58 ± 0.02, (14) 

GSC: vg — 0.62 ± 0.01. (15) 

Note that vq here practically coincides with the inverse fractal dimension of an open repulsive chain. Thus, we 
can interpret vq as the Flory swelling exponent, the value of which for large N is more accurately known from the 
renormalisation group calculation, namely vp — 0.588 ± 0.001. As discussed in Ref. [32] the GSC theory overestimates 
the Flory swelling exponent by about 0.04 as compared to the MC data for this range of finite polymer sizes. This is 
consistent with the amount by which the GSC values for vg exceed those of the MC simulation. 

At the same time, the prefactor ac does not obey a power law in G, but it can be fitted well via the following 
expression, 

GSC: a G = (0.74 ± 0.01) 2"^^-, MC: a G = (0.9 ± 0.03) 2~^n . (16) 
Or by using Eq. (^j again, we can express this for the good solvent via, 

3^ 1/4 



3ft 2 , ~ N uf „ £,1/4 ^-1/4 



-9 \ N 

ZR? g ~D v °(F-l)» D ( G+1 \ VD ^ VG -\. (18) 



Indeed, as can be seen from Tab. Ill, the resulting exponent vp> ~ 0.37 ± 0.01 in the GSC method, whereas from MC 
the resulting exponent value is vp> ~ 0.335 ± 0.008 > 1/3. These notations for the exponents vg, vp> at fixed values 
of either G or D are akin to those traditionally used in thermodynamics. 

For the poor solvent a similar procedure gives the following values for N > 200, 

GSC: Ag = 0.560 ± 0.007, v G = 0.331 ± 0.002 for [/(°) = 7, (19) 
GSC: Ag = 0.516 ± 0.004, v G = 0.333 ± 0.001 for [/(°) = 10, (20) 
MC: A G = 0.60 ±0.004, v G = 0.31 ± 0.02 for = 6, (21) 

or, roughly speaking, vq ~ v — |, with coefficient Ag decreasing slowly with increasing and being almost 
independent of other parameters such as G. Hence, we have for a poor solvent, 

3^2 „ N v _ D v a (p_ 1 j« (Gf+l) i VD=VG ^ V = }_ (22) 

It is interesting to note that for the ideal chain (absence of the excluded volume interactions, = 0) we obtain 
the following accurate asymptotic expressiorEj for large D and G, 



ZK 2 ~iVZD 1 / 2 G 1/2 . (23) 



Therefore, generally, the dependence of the radius of gyration on the spacer exponent is always the same as for 
an open/ring chain, but the ideal system lacks the complicated dependence on the number of generations G. It is 
worthwhile to comment also that such ideal system should not be confused with the system in the theta-solvent near 
the transition point J7g? . For the latter we shall have a scaling akin to Eqs. (|l8]]2|), 

3^ „ D *c {F _ 1) v D (G+i) ) Vg= 1 5 VJ} < VD < VD . (24) 



The values for v^, obtained from the GSC theory are reported in Tabs. |nj,|nj. We can conclude that indeed 
= 1/2 with a high accuracy, whereas v^> ~ 0.35. We expect the latter to be somewhat lower in the MC simulation 
than in the GSC theory again. 
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V. DISCUSSION AND CONCLUSION 



In this paper we have studied a range of homo- dendrimers with D = 1 — 6 spacers and G = 1 — 7 generations in 
the good and poor solvents, as well as across the coil-to-globule transition. These studies have been performed by 
two different techniques yielding a satisfactory agreement between their results. The first technique is based on the 
Gaussian self-consistent (GSC) method, while the second one is that of Monte Carlo (MC) simulation in continuous 
space. Both rely on the same coarse-grained model involving harmonic springs for connected monomers and the 
Lennard-Jones pair-wise potential with a hard core part. The overall accord between these two rather distinct 
approaches is quite important due to the limitations of both techniques. The MC technique, while in principle exactly 
modelling a given system, in practice is always limited by the statistical accuracy and quality of equilibration for very 
large dendrimers computationally affordable. The GSC method, on the other hand, is about 10 3 faster for the studied 
range of dendrimcr sizes and it directly produces the ensemble averaged equilibrium observables. However, it is only 
an approximate technique, which has some intrinsic inaccuracies, for instance somewhat overestimating the size of 
polymers in the good solvent. 

We have proposed a new scaling law for the radius of gyration of dendrimers which involves the two exponents vq 
(at fixed G) and Vjj (at fixed D), 

K g ~ D ua (F - 1)^( G +D ~ jya-vn N »d^ (25) 

where vq = v, with v 3/5, 1/2, and 1/3 being the inverse fractal dimension of an open/ring chain in the good, 
ideal and poor solvents respectively. While this agrees with the conclusions of Refs. |2^j2^ for vq, we have obtained 
a rather different dependence on the number of generations. Namely, we have concluded that Vn ~ v — 1/4 ~ 0.338 



for the good solvent, and s» v = 1/3 for the globule. This is consistent with Ref. 27, but our estimate of vjj for 
the good solvent is more accurate enabling us to claim that vjj is certainly above, although quite close to the poor 
solvent value 1/3. 

It is important to emphasise also that the law in Eq. ( p5j ) is valid for the model with the harmonic bonds employed 
here. Dendrimers with bonds possessing a finite stretchability and angular potentials would in all likelihood have 
more complicated and non-universal features in their behaviour. However, the harmonic model is an important and 
well tractable reference point for further studies at the level of full atomistic detail for concrete dendrimers, which 
could have bonds stretchable in different ranges. 

We have concluded that the coil-to-globule transition for dendrimers is continuous and that for the theta-solvent the 
scaling law for lZ g differs from that of the ideal dendrimer, in which lZ g ~ (DG) 1 ^ 2 , having instead the G-dependence 
of the form given by Eq. @ with 1/3 < v D < 3/5 - 1/4. 

Details of the conformational structure of dendrimers have been analysed for the good and poor solvents also. 
In some respects, our results support certain observations of Ref. ^4] made for the good solvent. Namely, we find 
a strong evidence of the re-entry of the terminal generation into the dendrimer up to its very centre, producing a 
rather delocalised partial density for those monomers. This results in the highest total density near the core for 
relatively small dendrimers, which decreases in larger dendrimers due to the 'entropic pull' effect of extra branches. 
However, as the number of generations exceeds G = 5 there also develops a lower density domain inside the dendrimer. 
Moreover, we have analysed the bond stretching along a continuous branch. This is highest near the core, especially 
for dendrimers with a large number of generations. The stretching decreases with increasing generation index and 
oscillates in the magnitude with the spacer index. This observation is interpreted consistently with the connectivity 
contribution to the mean energy. 

Furthermore, we have attempted to elucidate the degree of the steric congestion in a dendrimer by considering the 
shell index dependence of the monomer packing coefficient r\ p . The latter is the volume ratio of the space occupied 
by all other monomers around a given one. In a relatively small dendrimer r\ p oscillates in the spacer index with only 
a weak generation index dependence, except for the terminal generation which is less congested. A dendrimer with a 
large number of generations, however, has a less congested hollow domain inside it and the highest congestion in the 
penultimate generations in accord with the density distributions. 

The degree of asphericity of dendrimers appears to decrease with the increasing number of generations, depending 
rather weakly on the number of spacers. This is also confirmed by the behaviour of the static structure factor, which 
develops an oscillating globule-like character for larger G. 



Interestingly, while incorrect in concrete detail, the onion model predictions of Ref. 13 do have a rather limited 
meaning in a sense that there exists a hollow domain inside a dendrimer of large enough number of generations, but 
not at the core itself. However, the partial densities of generations fall off only beyond the appropriate cut-off radii. 
These densities are totally delocalised until the cut-offs (see Fig ||b), being fairly large at the core£3, and certainly 
are not concentrated within thin layers as surmised. 
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A similar analysis has been carried out for the poor solvent as well. The globular conformation of a dendrimer is 
much more compact, and it is rather uniform without any hollow domains. The packing coefficient for the globule is 
much larger in value than for the good solvent and depends on the spacer index only for the monomers in the inner 
generations. There is also much less stretching of bonds involved near the core than for the good solvent, even in a 
dendrimer with a large G. Generally, the bond stretching in the globule is stronger for the branch points than for the 
middle of spacers or terminal monomers. Thus, the connectivity of a dendrimer in the globular state manifests itself 
in a much more localised manner, and with a significantly smaller effect than for the good solvent. 

We believe that our results allow one to reconcile some of the disparate predictions from previous theoretical and 
simulation studies by finding some simple limiting regimes and by discovering non-monotonic dependencies of various 
observables and the local structure of dendrimers on the number of generations. These unusual features of dendrimers 
indeed make them extremely appealing for numerous future applications ranging from catalysis, rheology modifiers, 
to drug delivery and biotechnology more generally. 
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APPENDIX A: THE TOPOLOGICAL TREE 



The concept of the kinematic symmetries can be best explained by using the topological tree presented in Fig. 
|l6| , which is equivalent to the dendrimer diagram of Fig. [lj Hence, the horizontal line would represent any selected 
consecutive branch of the dendrimer, which we shall name the primary branch B, with the core monomer being on 
the left. The length of B is then ng = D(G + 1), excluding the core monomer itself. Any two branches become 
topologically non-equivalent to each other as long as they diverge from the common parent sub-branch at some value 
of the generation index g. We shall call any two such branches non-equivalent of order g. The leftmost vertical 
line on the left of Fig. fiq represents any branch non-equivalent to B of order g = 0, which we shall denote as Bo 
with the length ng = tiq. Similarly then, we shall have branches B g non-equivalent to B of order g with the length 
n g = D(G — p + l) for all values of g = 0, 1, ... G in our topological tree. Any two pairs of monomers then would have 
T>ij = T>iiji as long as the pair lies in the same places in the topological tree as the pair thereby defining 

the equivalence classes of monomer pairs. 

Moreover, we can easily count the number of independent degrees of freedom in the matrix 2?y. For that we 
have to count the number of possible pairs (p, p') in Fig. [l6] such that e.g. p' < p, where i — (p, ip) and i' — (p' , tp'). 
Accidentally, we note that the topological tree disregards the angular coordinates ip of any equivalent monomers. There 
are two possibilities of pairing p and p ': a) p,p' 6 B and b) p € B and p' e B g for some g. In the former case we get 
the number of pairs Cg = J27=i * = n B( n B + l)/2, whereas in the latter case we obtain C g — Yli= n * = n s( n s + 
Thus, we have derived the total number of independent degrees of freedom, 



g=G 



which is identical to Eq. (A-l) in Ref. 19 



Similarly, the total number of independent 1-point observables (e.g. the packing coefficient r]i), which depend on p 
only, is simply the length of the primary branch Bq, 

c\n d = D(G + 1) + 1. (A2) 
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APPENDIX: FIGURE CAPTIONS 



FIG. 1. Schematic connectivity diagram of the dendrimer D2G3. Here dashed circles denote different generations, big filled 
circles correspond to the branch points and the triangle is the core monomer. Notations for the indices d, p, <p and for some 
selected monomers, which are numbered in the clock-wise and outwards manner, are also introduced. 

FIG. 2. Plots of the mean-squared radius of gyration 31Z^ vs the number of generations G for the dendrimers in the good 
solvent (U^ ' = 0) from the GSC theory (filled circles and thin lines) and from the MC simulation (empty circles and thick 
lines). The curves (from bottom to top) correspond to the following values: D = 1 (GSC), 2 (MC), 2 (GSC), 4 (MC), 4 (GSC), 
6 (GSC). The largest relative statistical error, SR 2 g = AR 2 g /R 2 g ■ 100%, in this plot occurs for the D2G7 dendrimer and equals 
0.1%. 

FIG. 3. Plots of the mean ratios of the eigenvalues of the shape tensor A'™' vs the number of generations G for the 
dendrimers in the good solvent (C/' ' = 0) from MC simulation. Solid curves correspond to D = 2 and empty circles to D = 4. 
The largest relative statistical errors, 8\ (a) = AA (a) /A (a) • 100%, in this plot occur for the D2G7 dendrimer and are equal to: 
<5A« = 0.15%, = 0.2% and = 0.3%. 



FIG. 4. Typical snapshots of conformations of the dendrimers D3G3 (Fig. a) and D3G7 (Fig. b) in the good solvent from 
MC simulation. The core monomer is represented as a sphere of the monomer real diameter = I in both figures. 



FIG. 5. Kratky plots of the rescaled static structure factor S = p 2 31Z 2 S(p)/N vs the rescaled wave number p = p y^STZ 2 
for dendrimers in the good solvent (U^ = 0) from MC simulation. Here the solid line corresponds to D1G5, the long-dashed 
curve — to D1G6 and the short-dashed curve — to D2G5 dendrimer. The largest relative statistical errors in this plot occur 
for the D2G5 dendrimer and are equal to: Sp = 0.1% and 8S < 0.3%. 

FIG. 6. Plots of the monomer densities g g ^ vs the separation from the centre-of-mass r cm for the dendrimers D2G4 
(Fig. a) and D2G6 (Fig. b) in the good solvent (U^ = 0) from MC simulation. Insets of figures show the monomer densities 
of the same dendrimers in the poor solvent (U^ = 6). Curves from top to bottom as on the left-hand side correspond to 
the following densities. In main bodies of figures: the total density g g ^ (solid line), density of non-terminal generations 
^2 g<G 9(p (dashed line), density of the innermost generations + g^ (thin solid line with quadrangles), density of the 
terminal generation gQ (dotted line), and density of the penultimate generation g G V) _ 1 (thin line with triangles); In insets of 
figures: the total density g g ^ (solid line), density of non-terminal generations ^2 g<G g^ (dashed line), density of the 
terminal generation g G (dotted line) . Errorbars represent statistical errors of the data from MC simulation. 

FIG. 7. Plots of the mean-squared distances -D( p -i,o) O,o) between radial nearest neighbours vs the shell index p from MC 
simulation (main body of the figure) and from the GSC theory (inset of the figure) for dendrimers with D = 2 spacers in the 
good solvent (U (0) = 0). The curves correspond to G = 1,2,3,4,5,6,7 (from bottom to top). Errorbars represent statistical 
errors of the data from MC simulation. 

FIG. 8. Plots of the mean-squared distances Dt p _ 1 m ( p m vs the shell index p from MC simulation (main body of the figure) 
and from the GSC theory (inset of the figure) for dendrimers with D — 4 spacers in the good solvent (U (0) = 0). The curves 
correspond to G = 1, 2, 3, 4, 5, 6 (from bottom to top). Errorbars represent statistical errors of the data from MC simulation. 

FIG. 9. Plots of the mean-squared distances -D( p -i,o) ( p ,o) vs the shell index p from the GSC theory for dendrimers with 
D = 6 spacers in the good solvent (I7 (0) = 0). The curves correspond to G = 1, 2, 3, 4, 5, 6 (from bottom to top). 
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FIG. 10. Plots of the packing coefficient r\ p defined by Eqs. vs the rescaled shell index p/D for dendrimers consisting 

of G = 6 generations in the good solvent ([A ' = 0) from the GSC theory. The curves correspond to the following values of 
D — 1, 2, 3, 4, 6 (from top to bottom). 

FIG. 11. Plots of the packing coefficient rj p vs the shell index p for dendrimers with D = 3 spacers in the good solvent 
(U m = 0) from the GSC theory. The curves correspond to the following values of G: G = 1 (dotted line), G = 3 (short-dashed 
line), G = 5 (long-dashed line), and G — 7 (solid line). 

FIG. 12. Plots of the mean-squared radius of gyration 31Z 2 (Fig. a) and the specific energy slope N~ 1 d£/dU^ (Fig. b) 
vs the degree of Lennard-Jones attraction [/' ' across the coil-to-globule transition from the GSC theory for dendrimers with 
G — 5 generations and D = 1, 2, 3, 4, 6 (from bottom to top as on the left-hand side). 

FIG. 13. Plots of the mean-squared radius of gyration 2>TZ 2 g vs the number of generations G for the dendrimers in the poor 
solvent, = 7, from the GSC theory (filled circles and solid lines) and for the ideal dendrimers, d' ' = (filled quadrangles 
and dashed lines). The curves correspond to the following values of D for each of the two families: D — 1,2,3,4,5,6 (from 
bottom to top). 

FIG. 14. Plots of the mean-squared distances _D( p _ 10 ) ( p ,o) vs the shell index p from the GSC theory for dendrimers with 
D = 6 spacers in the poor solvent (U^ = 7). The curves correspond to G = 1, 2, 3, 4, 5, 6 (from bottom to top). 

FIG. 15. Plots of the packing coefficient r\ p vs the shell index p for dendrimers with D = 6 spacers in the poor solvent 
(j/(°> = 7) f ro m the GSC theory. The curves correspond to the following values of G = 1, 2, 3, 4, 5, 6 (from bottom to top). 

FIG. 16. Tree diagram topologically equivalent to that of Fig. [l] for D2G3 dendrimer. Dotted line illustrates the area 
satisfying condition p' < p. 
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APPENDIX: TABLES 



TABLE I. Values of the pair-wise interaction parameter at the coil-to-globule transition point Ug, and at the special point 
in parentheses vs D and G. The former has been determined by a quadratic interpolation of the maxima positions of 
—d 2 £/dU^ 2 . The latter has been determined by the condition, 1Z 2 (U^) = lZ 2 g (ideal). 





G= 1 


G = 2 


G = 3 


G = 4 


G= 5 


G = 6 


D = 


1 


4.86 (4.28) 


3.87 (3.81) 


3.27 (3.75) 


2.97 (4.02) 


2.75 (4.78) 


2.54 (6.97) 


D = 


2 


4.05 (3.50) 


3.24 (3.20) 


2.87 (3.17) 


2.70 (3.33) 


2.61 (3.77) 


2.50 (4.90) 


D = 


3 


3.56 (3.17) 


2.96 (2.95) 


2.68 (2.92) 


2.55 (3.03) 


2.48 (3.36) 


2.43 (4.14) 


D = 


4 


3.29 (2.99) 


2.80 (2.79) 


2.57 (2.77) 


2.45 (2.87) 


2.39 (3.13) 


2.36 (3.74) 


D = 


5 


3.11 (2.86) 


2.69 (2.70) 


2.49 (2.67) 


2.39 (2.75) 


2.35 (2.98) 




D = 


6 


2.98 (2.77) 


2.61 (2.62) 


2.43 (2.60) 


2.35 (2.68) 


2.31 (2.87) 





TABLE II. Values of parameters aa, vg and ac vg obtained from fitting procedure via Eqs. (|ll|Jl2J) for the good athermal, 
= 0, and the theta, , solvents over various values of D at a given value of G from the GSC and MC methods. 



G 


GSC: a G 


VG 


MC: a G 




GSC: a G 


VG 


1 


0.532 ±0.009 


0.615 ±0.005 


0.63 ± 0.01 


0.55 ±0.01 


0.60 ± 0.02 


0.477 ± 0.009 


2 


0.445 ± 0.007 


0.626 ± 0.004 


0.50 ±0.03 


0.58 ± 0.02 


0.51 ±0.02 


0.489 ± 0.006 


3 


0.371 ± 0.005 


0.629 ± 0.003 


0.44 ± 0.03 


0.57 ±0.01 


0.44 ±0.02 


0.500 ± 0.005 


4 


0.312 ±0.004 


0.628 ± 0.002 


0.35 ±0.01 


0.58 ±0.01 


0.39 ±0.01 


0.498 ± 0.004 


5 


0.266 ± 0.004 


0.624 ± 0.002 


0.29 ±0.02 


0.59 ± 0.01 


0.34 ±0.01 


0.501 ±0.006 


6 


0.232 ±0.003 


0.621 ±0.002 


0.26 ±0.01 


0.57 ±0.01 


0.30 ±0.01 


0.497 ± 0.005 



TABLE III. Values of parameters ao, vd and an, vd obtained from fitting procedure via Eqs. ( |ll| , ft2[ ) for the good athermal, 
— 0, and the theta, U^V , solvents over various values of G at a given value of D from the GSC and MC methods. 



D 


GSC: a D 


vd 


MC: a D 


v D 


GSC: a D 




1 


0.98 ±0.01 


0.377 ± 0.002 


1.08 ±0.01 


0.333 ± 0.002 


0.69 ± 0.03 


0.380 ± 0.006 


2 


1.18 ± 0.02 


0.373 ± 0.002 


1.29 ±0.02 


0.329 ± 0.002 


0.85 ±0.03 


0.355 ± 0.005 


3 


1.33 ±0.03 


0.370 ± 0.003 


1.43 ±0.03 


0.330 ± 0.003 


0.92 ±0.01 


0.346 ± 0.003 


4 


1.42 ±0.05 


0.370 ± 0.004 


1.50 ±0.05 


0.332 ± 0.004 


0.96 ±0.02 


0.346 ± 0.003 


5 


1.52 ±0.06 


0.369 ± 0.005 


1.48 ±0.09 


0.340 ± 0.007 


0.97 ±0.04 


0.349 ± 0.006 


6 


1.60 ± 0.07 


0.368 ± 0.005 


1.62 ±0.09 


0.335 ± 0.008 


1.00 ±0.04 


0.351 ± 0.005 
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